Double Spin Asymmetry of Electrons from Heavy Flavor Decays 
in p+p Collisions at y/s = 200 GeV 



j.o. mil, it. a. Mollis, vv. Moizmann, i\. Momma, a. Mong, i. Moragucm. i. Mori, u. MornDactt, 
3. Huang, 61 T. Ichihara, 50 * 51 R. Ichimiya, 50 H. Iinuma, 28 Y. Ikeda, 50 * 52 ' 60 K. Imai, 31 * 50 M. Inaba, 60 A. Iordanova, 6 

D. Isenhower, 1 M. Ishihara, 50 M. Issah, 61 A. Isupov, 26 D. Ivanischev, 49 Y. Iwanaga, 20 B.V. Jacak, 57 'Q J. Jia, 5 ' 56 
X. Jiang, 34 D. John, 58 B.M. Johnson, 5 T. Jones, 1 K.S. Joo, 41 D. Jouan, 47 J. Kamin, 57 S. Kaneti, 57 B.H. Kang, 19 

J.H. Kang, 64 J.S. Kang, 19 J. Kapustinsky, 34 K. Karatsu, 31 ' 50 M. Kasai, 50 ' 52 D. Kawall, 38 ' 51 A.V. Kazantsev, 30 
T. Kempel, 25 A. Khanzadeev, 49 K.M. Kijima, 20 B.I. Kim, 29 D.J. Kim, 27 E.-J. Kim, 8 Y.-J. Kim, 22 Y.K. Kim, 19 
E. Kinney, 11 A. Kiss, 15 E. Kistenev, 5 D. Kleinjan, 6 P. Kline, 57 L. Kochenda, 49 B. Komkov, 49 M. Konno, 60 

J. Koster, 22 D. Kotov, 49 A. Krai, 13 G.J. Kunde, 34 K. Kurita, 50 ' 52 M. Kurosawa, 50 Y. Kwon, 64 G.S. Kyle, 44 
R. Lacey, 56 Y.S. Lai, 12 J.G. Lajoie, 25 A. Lebedev, 25 D.M. Lee, 34 J. Lee, 16 K.B. Lee, 29 K.S. Lee, 29 S.H. Lee, 57 

S.R. Lee, 8 M.J. Leitch, 34 M.A.L. Leite, 54 X. Li, 9 S.H. Lim, 64 L.A. Linden Levy, 11 A. Litvinenko, 26 H. Liu, 34 
MX. Liu, 34 B. Love, 61 D. Lynch, 5 C.F. Maguire, 61 Y.I. Makdisi, 4 A. Malakhov, 26 A. Manion, 57 V.I. Manko, 30 

E. Mannel, 12 Y. Mao, 48 ' 50 H. Masui, 60 M. McCumber, 57 P.L. McGaughey, 34 D. McGlinchey, 17 C. McKinney, 22 
N. Means, 57 M. Mendoza, 6 B. Meredith, 22 Y. Miake, 60 T. Mibe, 28 A.C. Mignerey, 37 K. Miki, 50 > 60 A. Milov, 62 

J.T. Mitchell, 5 Y. Miyachi, 50 * 59 A.K. Mohanty, 3 H.J. Moon, 41 Y. Morino, 10 A. Morreale, 6 D.P. Morrison, 5 
S. Motschwiller, 40 T.V. Moukhanova, 30 T. Murakami, 31 J. Murata, 50 ' 52 S. Nagamiya, 28 J.L. Nagle, 11 M. Naglis, 62 
M.I. Nagy, 63 I. Nakagawa, 50 ' 51 Y. Nakamiya, 20 K.R. Nakamura, 31 ' 50 T. Nakamura, 50 K. Nakano, 50 J. Newby, 33 
M. Nguyen, 57 M. Nihashi, 20 R. Nouicer, 5 A.S. Nyanin, 30 C. Oakley, 18 E. O'Brien, 5 CA. Ogilvie, 25 M. Oka, 60 
K. Okada, 51 A. Oskarsson, 36 M. Ouchida, 20 ' 50 K. Ozawa, 10 R. Pak, 5 V. Pantuev, 23 ' 57 V. Papavassiliou, 44 
B.H. Park, 19 I.H. Park, 16 S.K. Park, 29 S.F. Pate, 44 H. Pei, 25 J.-C. Peng, 22 H. Pereira, 14 V. Peresedov, 26 
DYu. Peressounko, 30 R. Petti, 57 C. Pinkenburg, 5 R.P. Pisani, 5 M. Proissl, 57 M.L. Purschke, 5 H. Qu, 18 J. Rak, 27 
I. Ravinovich, 62 K.F. Read, 46 ' 58 K. Reygers, 39 V. Riabov, 49 Y. Riabov, 49 E. Richardson, 37 D. Roach, 61 
G. Roche, 35 S.D. Rolnick, 6 M. Rosati, 25 S.S.E. Rosendahl, 36 P. Rukoyatkin, 26 B. Sahlmueller, 39 ' 57 N. Saito, 28 
T. Sakaguchi, 5 V. Samsonov, 49 S. Sano, 10 M. Sarsour, 18 T. Sato, 60 M. Savastio, 57 S. Sawada, 28 K. Sedgwick, 6 
R. Seidl, 51 R. Seto, 6 D. Sharma, 62 I. Shein, 21 T.-A. Shibata, 50 * 59 K. Shigaki, 20 H.H. Shim, 29 M. Shimomura, 60 
K. Shoji, 31 ' 50 P. Shukla, 3 A. Sickles, 5 C.L. Silva, 25 D. Silvermyr, 46 C. Silvestre, 14 K.S. Sim, 29 B.K. Singh, 2 
CP. Singh, 2 V. Singh, 2 M. Slunecka, 7 T. Sodre, 40 RA. Soltz, 33 W.E. Sondheim, 34 S.P. Sorensen, 58 
I.V. Sourikova, 5 P.W. Stankus, 46 E. Stcnlund, 36 S.P. Stoll, 5 T. Sugitate, 20 A. Sukhanov, 5 J. Sun, 57 J. Sziklai, 63 
E.M. Takagui, 54 A. Takahara, 10 A. Taketani, 50 * 51 R. Tanabe, 60 Y. Tanaka, 42 S. Taneja, 57 K. Tanida, 55 ' 31 > 50 
M.J. Tannenbaum, 5 S. Tarafdar, 2 A. Taranenko, 56 E. Tennant, 44 H. Thcmann, 57 D. Thomas, 1 M. Togawa, 51 

L. Tomasek, 24 M. Tomasek, 24 H. Torii, 20 R.S. Towell, 1 I. Tserruya, 62 Y. Tsuchimoto, 20 K. Utsunomiya, 10 
C. Vale, 5 H.W. van Hecke, 34 E. Vazquez-Zambrano, 12 A. Veicht, 12 J. Velkovska, 61 R. Vertesi, 63 M. Virius, 13 
A. Vossen, 22 V. Vrba, 24 E. Vznuzdaev, 49 X.R. Wang, 44 D. Watanabe, 20 K. Watanabe, 60 Y. Watanabe, 50 ' 51 
Y.S. Watanabe, 10 F. Wei, 25 R. Wei, 56 J. Wessels, 39 S.N. White, 5 D. Winter, 12 C.L. Woody, 5 R.M. Wright, 1 
M. Wysocki, 11 Y.L. Yamaguchi, 10 R. Yang, 22 A. Yanovich, 21 J. Ying, 18 S. Yokkaichi, 50 ' 51 J.S. Yoo, 16 Z. You, 34 ' 48 



2 



G.R. Young, 46 I. Younus, 43 I.E. Yushmanov, 30 W.A. Zajc, 12 A. Zelenski 4 S. Zhou, 9 and L. Zolin 26 

(PHENIX Collaboration) 

Abilene Christian University, Abilene, Texas 79699, USA 
2 Department of Physics, Banaras Hindu University, Varanasi 221005, India 
s Bhabha Atomic Research Centre, Bombay 400 085, India 
4 Collider- Accelerator Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
5 Physics Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
6 University of California - Riverside, Riverside, California 9252,1, USA 
'Charles University, Ovocny trh 5, Praha 1, 116 36, Prague, Czech Republic 
8 Chonbuk National University, Jeonju, 561-756, Korea 
9 Science and Technology on Nuclear Data Laboratory, China Institute of Atomic Energy, Beijing 1024-13, P. R. China 
10 Center for Nuclear Study, Graduate School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

11 University of Colorado, Boulder, Colorado 80309, USA 
12 Columbia University, New York, New York 10027 and Nevis Laboratories, Irvington, New York 10533, USA 
13 Czech Technical University, Zikova 4, 166 36 Prague 6, Czech Republic 
14 Dapnia, CEA Saclay, F-91191, Gif-sur-Yvette, France 
15 ELTE, Eotvos Lordnd University, H - 1117 Budapest, Pdzmdny P. s. 1/A, Hungary 
16 Ewha Womans University, Seoul 120-750, Korea 
"Florida State University, Tallahassee, Florida 32306, USA 
18 Georgia State University, Atlanta, Georgia 30303, USA 
19 Hanyang University, Seoul 133-792, Korea 
20 Hiroshima University, Kagamiyama, Higashi-Hiroshima 739-8526, Japan 
21 IHEP Protvino, State Research Center of Russian Federation, Institute for High Energy Physics, Protvino, 142281, Russia 

22 University of Illinois at Urbana- Champaign, Urbana, Illinois 61801, USA 
23 Institute for Nuclear Research of the Russian Academy of Sciences, prospekt 60-letiya Oktyabrya 7a, Moscow 117312, Russia 
24 Institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague 8, Czech Republic 

25 Iowa State University, Ames, Iowa 50011, USA 
26 Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 
27 Helsinki Institute of Physics and University of Jyvaskyld, P.O.Box 35, FI-40014 Jyvdskyld, Finland 
28 KEK, High Energy Accelerator Research Organization, Tsukuba, Ibaraki 305-0801, Japan 
29 Korea University, Seoul, 136-701, Korea 
30 Russian Research Center "Kurchatov Institute", Moscow, 123098 Russia 
31 Kyoto University, Kyoto 606-8502, Japan 
32 Labor atoire Leprince-Ringuet, Ecole Polytechnique, CNRS-IN2P3, Route de Saclay, F-91128, Palaiseau, France 
33 Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
34 Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
35 LPC, Universite Blaise Pascal, CNRS-IN2P3, Clermont-Fd, 63177 Aubiere Cedex, France 
36 Department of Physics, Lund University, Box 118, SE-221 00 Lund, Sweden 
37 University of Maryland, College Park, Maryland 20742, USA 
38 Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003-9337, USA 
39 Institut fur Kernphysik, University of Muenster, D-48149 Muenster, Germany 
40 Muhlenberg College, Allentown, Pennsylvania 18104-5586, USA 
41 Myongji University, Yongin, Kyonggido 449-728, Korea 
42 Nagasaki Institute of Applied Science, Nagasaki-shi, Nagasaki 851-0193, Japan 
43 University of New Mexico, Albuquerque, New Mexico 87131, USA 
44 New Mexico State University, Las Cruces, New Mexico 88003, USA 
45 Department of Physics and Astronomy, Ohio University, Athens, Ohio 45701, USA 
46 Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
47 IPN-Orsay, Universite Paris Sud, CNRS-IN2P3, BP1, F-91406, Orsay, France 
48 Peking University, Beijing 100871, P. R. China 
49 PNPI, Petersburg Nuclear Physics Institute, Gatchina, Leningrad region, 188300, Russia 
5 "RIKEN Nishina Center for Accelerator-Based Science, Wako, Saitama 351-0198, Japan 
51 RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
52 Physics Department, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima, Tokyo 171-8501, Japan 
53 Saint Petersburg State Polytechnic University, St. Petersburg, 195251 Russia 
54 Universidade de Sao Paulo, Instituto de Fisica, Caixa Postal 66318, Sao Paulo CEP05315-970, Brazil 

J Seoul National University, Seoul, Korea 
56 Chemistry Department, Stony Brook University, SUNY, Stony Brook, New York 11794-3400, USA 
57 Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, New York 11794-3400, USA 

58 University of Tennessee, Knoxville, Tennessee 37996, USA 
59 Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152-8551, Japan 
60 Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
61 Vanderbilt University, Nashville, Tennessee 37235, USA 



3 



Weizmann Institute, Rehovot 76100, Israel 
63 Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Hungarian 
Academy of Sciences (Wigner RCP, RMKI) H-1525 Budapest 114, POBox 49, Budapest, Hungary 
64 Yonsei University, IPAP, Seoul 120-749, Korea 
(Dated: September 17, 2012) 

We report on the first measurement of double-spin asymmetry, All, of electrons from the decays 
of hadrons containing heavy flavor in longitudinally polarized p+p collisions at \fs = 200 GeV for 
Pt = 0.5 to 3.0 GeV /c The asymmetry was measured at midrapidity (|?j| < 0.35) with the PHENIX 
detector at the Relativistic Heavy Ion Collider. The measured asymmetries are consistent with zero 
within the statistical errors. We obtained a constraint for the polarized gluon distribution in the 
proton of \Ag/g(log 10 x — —1.6t.o'\, H = rn^)\ 2 < 0.033 (la) based on a leading-order perturbativc- 
quantum-chromodynamics model, using the measured asymmetry. 

PACS numbers: 13.85.Ni, 13.88. +e,14.20.Dh,25.75.Dw 



I. INTRODUCTION 

The measurement of the first moment of the proton's 
spin-dependent structure function g\ by the European 
Muon Collaboration (EMC) GIB revealed a discrepancy 
from the Ellis- Jaffe sum rule [3|, [j] and also the fact that 
the SU(3) flavor-singlet axial charge g^ 1 was smaller than 
expected from the static and relativistic quark models Q . 
After these discoveries, experimental efforts [H-HI focused 
on a detailed understanding of the spin structure of the 
proton. The proton spin s z /Ti — 1/2 can be decomposed 



as 



±A£ + AG 



L z from conservation of angular 
momentum. The measurements precisely determined the 
total spin carried by quarks and anti-quarks, AS, which 
is only about 30% of the proton spin. The remaining pro- 
ton spin can be attributed to the other components, the 
gluon spin contribution (AG) and/or orbital angular mo- 
mentum contributions (L z ). The total gluon polarization 
is given by 



(1) 



AG(^) = / dxAg(x,n), 



where x and /i represent Bjorken x and factorization scale 
respectively. The challenge for the AG(fi) determina- 
tion is to precisely map the gluon polarization density 
Ag(x, /i) over a wide range of x. 

The Relativistic Heavy Ion Collider (RHIC), which can 
accelerate polarized proton beams up to 255 GeV, is a 
unique and powerful facility to study the gluon polariza- 
tion. One of the main goals of the RHIC physics program 
is to determine the gluon polarization through measure- 
ments of longitudinal double-spin asymmetries, 



All = 



(2) 



where a ++ and er + denote the cross sections of a spe- 
cific process in the polarized p+p collisions with same 



and opposite helicities. Using All, the polarized cross 
sections, <r ++ and er H , can be represented as, 



<j +± =a {l±A LL ), 



(3) 



'Deceased 

^PHENIX Spokesperson: jacak@skipper.physics.sunysb.edu 



where is the unpolarized cross section of the process. 
All has been measured previously in several channels 
by PHENIX and STAR, including inclusive tt° [§-[13], -q 
[13], and jet [H4il] production. 

Using the measured asymmetries, as well as the world- 
data on polarized inclusive and semi-inclusive deep- 
inelastic scattering (DIS) [13, EH, a global anal- 
ysis based on perturbative-quantum-chromodynamics 
(pQCD) calculation was performed at next-to-leading or- 
der (NLO) in the strong-coupling constant as [l9| . The 
resulting Ag(x, /i) from the best fit is too small to explain 
the proton spin in the Bjorken x range of 0.05 < x < 0.2 
(—1.3 < log 10 x < —0.7) without considering L Z) though 
a substantial gluon polarization is not ruled out yet due 
to the uncertainties. Also, due to the limited Bjorken x 
coverage, there is a sizable uncertainty in Eq. [I] from the 
unexplored small x region. 

The polarized cross section of heavy flavor production 
on the partonic level is well studied with leading-order 
(LO) and NLO pQCD calculations [UH. The heavy 
quarks are produced dominant ly by the gluon-gluon in- 
teraction at the partonic level [23|. Therefore, this chan- 
nel has good sensitivity to the polarized gluon density. 
In addition, the large mass of the heavy quark ensures 
that pQCD techniques are applicable for calculations of 
the cross section. Therefore, the measurement of heavy 
flavor production in polarized proton collisions is a useful 
tool to study gluon polarization. 

In p+p collisions at y/s = 200 GeV, the heavy flavor 
production below p? ~ 5 GeV/c is dominated by charm 
quarks. The Bjorken x region covered by this process 
at midrapidity is centered around 2m c /y/s ~ 1.4 x 10~ 2 
where m c represents the charm quark mass. Hence, mea- 
surement of the spin dependent heavy flavor production 
is sensitive to the unexplored x region, and complements 
other data on the total gluon polarization AG(/i). 

At PHENIX, hadrons containing heavy flavors are 
measured through their semi-leptonic decays to electrons 
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and positrons (heavy flavor electrons) [24| HH . Therefore 
the double-spin asymmetry of the heavy flavor electrons 
is an important measurement for the gluon polarization 
study. In this paper, we report the first measurement of 
this asymmetry, and a resulting constraint on the gluon 
polarization with an LO pQCD calculation. 

The organization of this paper is as follows: We intro- 
duce the PHENIX detector system used for the measure- 
ment in Sec.HU The method for the heavy flavor electron 
analysis is discussed in Sec. lIIII and the results of the cross 
section and the spin asymmetry are shown in Sec. lIVI and 
Sec. El respectively. From the asymmetry result, we es- 
timate a constraint on the polarized gluon density, which 
is described in Sec. I VII For the sake of simplicity, we use 
the word "electron" to include both electron and positron 
throughout this paper, and distinguish by charge where 
necessary. 



II. EXPERIMENTAL SETUP 

This measurement is performed with the PHENIX de- 
tector positioned at one of collision points at RHIC. 
The RHIC accelerator comprises the blue ring circulat- 
ing clockwise and the yellow ring circulating counter- 
clockwise. For this experiment, polarized bunches are 
stored and accelerated up to 100 GeV in each ring and 
collide with longitudinal polarizations of ~ 57% along 
the beams at the collision point with a collision energy of 
y/s = 200 GeV. The bunch polarizations are changed to 
parallel (beam-helicity +) or anti-parallel (beam-helicity 
— ) along the beams alternately in the collisions to real- 
ize all 4 (— 2 x 2) combinations of the crossing beam- 
helicities. Each time the accelerator is filled, the pattern 
of beam helicities in the bunches is changed, in order 
to confirm the absence of a pattern dependence of the 
measured spin asymmetry. See Sec. [V] for details. 

A detailed description of the complete PHENIX detec- 
tor system can be found elsewhere 26 32] . The main 
detectors that are used in this analysis are beam-beam 
counters (BBC), zero degree calorimeters (ZDC), and two 
central arm spectrometers. The BBC provides the col- 
lision vertex information and the minimum bias (MB) 
trigger. The luminosity is determined by the number of 
MB triggers. Electrons are measured with the two cen- 
tral spectrometer arms which each cover a pseudorapidity 
range of \rj\ < 0.35 and azimuthal angle Ac/> — ir/2. 

Figure Q] shows the beam view of the 2009 PHENIX 
central arms configuration, which comprises the central 
magnet (CM), drift chamber (DC), and pad chamber 
(PC) [for charged particle tracking], the ring- imaging 
Cerenkov detector (RICH) and hadron blind detector 
(HBD) [33j, l34| [for electron identification], and the elec- 
tromagnetic calorimeter (EMCal) [for energy measure- 
ment] . Below we summarize the features of the detectors 
and the CM. 

The BBCs are two identical counters positioned at 
±1.44 m from the nominal interaction point along the 



PHENIX Central Arm 




West Beam View East 



FIG. 1: (color online) Beam view (at z = 0) of the PHENIX 
central arm detectors in 2009. See text for details. 



beam direction and cover pseudorapidity of 3.1 < < 
3.9. They measure the collision vertex along the beam 
axis by measuring the time difference between the two 
counters, and also provide the MB trigger defined by at 
least one hit on each side of the vertex. The position 
resolution for the vertex is ~ 2.0 cm in p+p collision. 

The ZDCs, which are located at ±18.0 m away from 
the nominal interaction point along the beam direc- 
tion, detect neutral particles near the beam axis (9 < 
2.5 mrad). Along with the BBCs, the trigger counts 
recorded by the ZDCs are used to determine the rela- 
tive luminosity between crossings with different beam- 
helicities combinations. The ZDCs also serve for mon- 
itoring the orientation of the beam polarization in the 
PHENIX interaction region through the experiment. 

The transverse momentum of a charged particle track 
is determined by its curvature in the mag netic field pro- 
vided by the PHENIX CM system [27]. The CM is 
energized by two pairs of concentric coils and provides 
an axial magnetic field parallel to the beam direction. 
During this measurement, the two coils of the CM were 
operated in the canceling ("H — ") configuration. This 
configuration is essential for the background rejection of 
the heavy flavor electron measurement with the HBD 
as described later. In this configuration, the field is 
almost canceled out around the beam axis in the ra- 
dial region < R < 50 cm, and has a peak value of 
~ 0.35 T around R ~ 100 cm. The total field integral is 
\fBxdl\ = 0.43 Tm. 

The DC and PC in the central arms measure charged 
particle trajectories in the azimuthal direction to deter- 
mine the transverse momentum (pt) of each particle. By 
combining the polar angle measured by the PC and the 
vertex information along the beam axis from the BBC 
with px, the total momentum p is determined. The DC 
is positioned between 202 cm and 246 cm in radial dis- 
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tance from the collision point for both the west and east 
arms and the PC is 247-252 cm. 

The RICH is a threshold-type gas Cerenkov counter 
and the primary detector used to identify electrons in 
PHENIX. It is located in the radial region of 2.5-4.1 m. 
The RICH has a Cerenkov threshold of 7 = 35, which 
corresponds to p = 20 MeV/c for electrons and p = 
4.9 GeV/c for charged pions. 

The EMCal comprises four rectangular sectors in each 
arm. The six sectors based on lead-scintillator calorime- 
try and the two (lowest sectors on the east arm) based on 
lead-glass calorimetry are positioned at radial distances 
from the collision point of ~ 5.1 m and ~ 5.4 m, respec- 
tively. 




(a)HBD from top view. 




frame 

(b)HBD exploded view. 

FIG. 2: (color online) (a) Top view of the HBD showing the 
location of the HBD in the central magnet, (b) Exploded 
view of one half HBD arm. CF4 gas is filled in the volume as 
the Cerenkov light radiator. 



A challenging issue for the heavy flavor electron mea- 
surement is to reject the dominant background of elec- 
tron pairs from 7 conversions and Dalitz decays of 7r° 
and 77 mesons, which are mediated by virtual photons. 
These electrons are called "photonic electrons" , while all 
the other electrons are called "nonphotonic electrons". 
Most nonphotonic electrons are from heavy flavor decays, 
however, electrons from K e s decays (K — > evix) and the 
dielectron decays of light vector mesons are also nonpho- 
tonic 24]. The HBD aims to considerably reduce the 
photonic electron pairs utilizing distinctive feature of the 
e + e~ pairs, namely their small opening angles. 

The HBD is a position-sensitive Cerenkov detector op- 
erated with pure CF4 g radiator. It covers pseudo- 
rapidity \r\\ < 0.45 and 2 x 3ir/4 in azimuth. The cover- 
age is larger than the acceptance of the other detectors in 
the central arm in order to detect photonic electron pairs 
with only one track reconstructed in the central arm and 
the other outside of the central arm acceptance. Fig- 
ure shows the top view and exploded view of the HBD. 
The HBD has a 50 cm long radiator directly coupled in a 
windowless configuration to a readout element consisting 
of a triple Gas Electron Multiplier (GEM) stack, with a 
Csl photocathode evaporated on the top surface of the 
GEM facing the collision point and pad readout at the 
exterior of the stack. The readout element in each HBD 
arm is divided into five sectors. The expected number of 
photoelectrons for an electron track is about 20, which 
is consistent with the measured number. Since the HBD 
is placed close to the collision point, the material thick- 
ness is small in order to minimize conversions. The total 
thickness to pass through the HBD is 0.024Xrj and the 
thickness before the GEM pads is 0.007X . 

The Cerenkov light generated by electrons is directly 
collected on a photosensitive cathode plane, forming an 
almost circular blob image. The readout pad plane com- 
prises hexagonal pads with an area of 6.2 cm 2 (hexagon 
side length a = 1.55 cm) which is comparable to, but 
smaller than, the blob size which has a maximum area of 
9.9 cm 2 . 

The HBD is located in a field free region that preserves 
the original direction of the e + e~ pair. The Cerenkov 
blobs created by electron pairs with a small opening an- 
gle overlap, and therefore generate a signal in the HBD 
with twice the amplitude of a single electron. Electrons 
originating from 7r° and r\ Dalitz decays and 7 conver- 
sions can largely be eliminated by rejecting tracks which 
correspond to large signals in the HBD. 



III. HEAVY FLAVOR ELECTRON ANALYSIS 

With the improved signal purity from the HBD, the 
double helicity asymmetry of the heavy flavor electrons 
was measured. In this section, we explain how the heavy 
flavor electron analysis and the purification of the heavy 
flavor electron sample using the HBD was performed. 
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A. Data Set 

The data used here were recorded by PHENIX during 
2009. The data set was selected by a level- 1 electron 
trigger in coincidence with the MB trigger. The electron 
trigger required a minimum energy deposit of 0.6 GeV in 
a 2 x 2 tile of towers in EMCal, Cerenkov light detection 
in the RICH, and acceptance matching of these two hits. 
After a vertex cut of |z vtx | < 20 cm and data quality cuts, 
an equivalent of 1.4 x 10 11 MB events, corresponding to 
6.1 pb _1 , sampled by the electron trigger were analyzed. 



B. Electron Selection 

Electrons are reconstructed using the detectors in the 
PHENIX central arm described above. Several useful 
variables for the electron selection which were used in 
the previous electron analysis in 2006 [25[ are also used in 
this analysis. In addition to the conventional parameters, 
we introduced a new value, q c \ ns , for the HBD analysis. 

hbdeharge: q c \ us Total charge of the associated HBD 
cluster calibrated in units of the number of photo- 
electrons (p.e.). 

The electron selection cuts (elD-cut) are: 

4.0cr matching between track and EMCal cluster 

# of hit tubes in RICH around track > 2 
3.5cr matching between track and HBD cluster 
shower profile cut on EMCal 

0.57 < E/p < 1.37 (0.5 GeV/c < p T < 1.0 GeV/c) 
0.60 < E/p < 1.32 (1.0 GeV/c < p T < 1.5 GeV/c) 
0.64 < E/p < 1.28 (1.5 GeV/c < p T < 5.0 GeV/c) 

# of hit pads in HBD cluster > 2 
(jcius > 8.0 p.e. 

(<7cius > 4.0 p.e. for one low-gain HBD sector) 

These cuts require hits in the HBD, RICH, and EM- 
Cal that are associated with projections of the track onto 
these detectors. The shower profile in the EMCal is re- 
quired to match the profile expected of an electromag- 
netic shower. For electrons, the energy deposit on EM- 
Cal, E, and the magnitude of the reconstructed momen- 
tum on DC and PC, p, should match due to their small 
mass. Therefore the ratio, E/p, was required to be close 
to 1. Since the energy resolution of the EMCal depends 
on the momentum of the electron, the cut boundaries 
were changed in different momentum range. Charged 
particles traversing the CF4 volume in the HBD produce 
also scintillation light, which has no directivity and cre- 
ates hits with small charge in random locations in the 
GEM pads. To remove HBD background hits by the scin- 
tillation light, a minimum charge and a minimum cluster 
size were required for the HBD hit clusters. During this 
measurement, the efficiency for the Cerenkov light in one 
HBD sector was low compared with other sectors. Hence 



we apply a different charge cut to that HBD sector for 
the electron selection. 

The E/p distribution for tracks selected with these cuts 
is shown in Fig. [3] The clear peak around E/p = 1 corre- 
sponds to electrons and the spread of events around the 
peak consists mainly of electrons from K e ^ decays and 
misidentified hadrons. As the figure shows, the fraction 
of these background tracks in the reconstructed electron 
sample after applying elD-Cut including the E/p cut was 
small. The fractions of the K e % decays and the misiden- 
tified hadrons are described in Sec. IIIIDI and Sec. MI El 

As mentioned in Sec. [Ill we remove the photonic elec- 
trons and purify the heavy flavor electrons on the basis 
of the associated HBD cluster charge. The nonphotonic 
electron cuts (npe-Cut) are: 

8.0 < g c i us < 28.0 p.e. 

(4.0 < g c i us < 17.0 p.e. for 1 low-gain HBD sector) 



.xia" 




FIG. 3: E/p distributions for 0.5 GeV/c < p T < 1.0 GeV/c 
reconstructed charged tracks with the elD-Cut other than the 
E/p cut. Criteria of the E/p cut for the momentum region 
are shown by dashed lines in the plot. 



Yield estimation of heavy flavor electrons with 
HBD 



We categorize the HBD hit clusters into three types ac- 
cording to the source of the cluster. A cluster created by 
a single blob of Cerenkov light from a nonphotonic elec- 
tron as shown in Fig. 4(a) is defined as a single cluster. 



On the other hand, a cluster created by merging blobs of 
Cerenkov light from a track pair of photonic electrons as 
shown in Fig. |4(b)| is defined as a merging cluster. How- 
ever, a portion of the photonic electrons which have a 
large enough opening angle such that the two cluster do 
not merge (typically ^,0.1 r ad) cr eates two separated 
single clusters as shown in Fig. 4(c) Therefore the single 



clusters are created by both of the nonphotonic electron 
and the photonic electron with a large opening angle. 
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We also define another type of cluster created by scin- 
tillation light, which we call a scintillation cluster. Scin- 
tillation hits which accidentally have large hit charges 
and have neighboring hits can constitute clusters. Pho- 
tonic electrons from 7 conversions after the HBD GEM 
pads do not create Cerenkov light in the HBD gas volume. 
Hence they basically do not have associated clusters in 
the HBD and they are rejected by the HBD hit require- 
ment in the elD-Cut. However, a portion of these are 
accidentally associated with scintillation clusters and sat- 
isfy the elD-Cut and so also survive in the reconstructed 
electron sample. 

In Sec. IIII C 1[ we estimated yields of these clusters 
from the distribution shape of the HBD cluster charge. 
We also estimated the small component of single clus- 
ters generated from photonic electrons which have the 
large opening angles as described in Sec. IIII C 21 Then 
we determined the nonphotonic electron yield. Subtract- 
ing additional background electrons from K e 3 decays and 
e + e~ decays of light vector mesons, we obtain the heavy 
flavor electron yield as described in Sec. IIII C 31 




(a) Non-photonic electrons. 




(b) Photonic electrons (merging (c) Photonic electrons 

cluster). (separated clusters). 



FIG. 4: (color online) Responses of the HBD for (a) non- 
photonic electrons and (b,c) photonic electrons, (b) Most of 
the photonic electron pair create merging clusters, (c) How- 
ever, the photonic electrons with large opening angles create 
separated clusters. 



1. Yield estimation of single clusters 

All clusters associated with the reconstructed electrons 
can be classified into the above three types. The yield of 
the electrons associated with the single clusters must be 
evaluated to estimate the yield of the heavy flavor elec- 
trons. The shapes of the q c i us distributions for the three 
cluster types are quite different since merging clusters 



have basically double the charge of single clusters and 
the charge of scintillation clusters is considerably smaller 
than the charge of the single cluster. Using the differ- 
ence in the shapes, we estimate yields of these clusters as 
follows. 

The probability distributions of g c ius for single and 
merging clusters were estimated by using low-mass 
unlike-sign electron pairs reconstructed with only the 
elD-Cut, which is dominated by photonic electron pairs. 
We defined the unlike-sign electron pairs whose two elec- 
trons were associated with two different HBD clusters 
as separated electron pairs and the pairs whose two elec- 
trons were associated to the same HBD cluster as merging 
electron pairs. The probability distribution of q c i us for 
the single clusters were estimated by the q c \ us distribution 
of the separated electron pairs and the probability distri- 
bution of g c i us for the merging clusters were estimated by 
the Qcius distribution of the merging electron pairs. The 
reconstruction of the electron pairs creates a small bias 
on the shapes of the Q c ius distributions. Corrections for 
this bias are estimated by simulation and applied to the 
distributions. The probability distributions are denoted 
as /c(<7cius) for the single clusters and /™( Qcius) for the 
merging clusters. The probability distribution of q c i U s for 
the scintillation clusters is also estimated by the distri- 
bution of the hadron tracks reconstructed by the DC/PC 
tracking and the RICH veto and denoted as /c C1 (q c 1us)- 

The variables used in the hbdeharge analysis are: 

/,? Probability distribution of q c i U s for the single clusters 
/™ Probability distribution of q c i U s for the merging clus- 
ters 

ysci p ro bability distribution of q c i U s for the scintillation 
clusters 

n s Number of single clusters after applying elD-Cut. 
n m Number of merging clusters after applying elD-Cut. 
n sc i Number of scintillation clusters after applying elD- 
Cut. 

n s Number of single clusters after applying elD-Cut and 
npe-Cut. 

h m Number of merging clusters after applying elD-Cut 
and npe-Cut. 

n sc j Number of scintillation clusters after applying elD- 
Cut and npe-Cut. 

The Qcius distribution of the reconstructed electrons 
found by applying elD-Cut is fitted with a superposition 
of the three probability distributions 

n s X /c( Qcius) + 

n m x /-(Qcius) + (4) 

«sci X / c SCi (<7clus), 

where n s , n m and n SC i are fitting parameters that rep- 
resent the numbers of the reconstructed electrons asso- 
ciating to single clusters, merging clusters and scintilla- 
tion clusters after applying elD-Cut respectively. The 
fraction of nonphotonic electrons and photonic electrons 
are different in different px region of the reconstructed 



8 



electron sample. Therefore the fitting was performed for 
each pt region and u s {pt), n m (px) and n sc i(pT) for each 
Pt region were determined. In the fitting, the distribu- 
tion functions, / c s (g c i us ), /^(gtlus), and / c SCI (g c i us ), are as- 
sumed to be pt independent because the velocity of elec- 
trons in pt region of interest is close enough to the speed 
of light in vacuum such that the yield of Cerenkov light 
from the electron is nearly independent of pt- We also 
compared the shapes of the distributions in different pt 
regions to confirm that the effect from the track curvature 
is small enough to be ignored even at pt ~ 0.5 GeV/c. 
On the other hand, / c s ((? c ius), /™(<?clus) and / c SCI (g c ius) for 
different HBD sectors vary slightly. Considering this dif- 
ference, the fitting is performed for each sector individu- 
ally. 

The <7dus distribution for the reconstructed electrons 
with transverse momentum pt ranging from 0.75 GeV/c 
to 1.00 GeV/c and the fitting result are shown in Fig. [5] 
for one HBD sector. The charge distribution of the recon- 
structed electrons is well reproduced by the superposition 
of the three individual components. 
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merged cluster : n m 
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FIG. 6: (color online). Yield spectra of HBD clusters after 
applying elD-Cut and npe-Cut estimated from the HBD clus- 
ter charge fitting. The plot shows the spectrum for the single 
clusters (n a , blue circle), the spectrum for the merging clus- 
ters (n m , red square) and the spectrum for the scintillation 
clusters (n S ci, green triangle). The error bars represent fitting 
uncertainties. 
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FIG. 5: (color online). Charge distribution of HBD clusters 
associated with reconstructed electrons with a transverse mo- 
mentum ranging from 0.75 GeV/c to 1.00 GeV/c (solid black 
line), and the charge distribution for each component, i.e., 
single clusters (n s f^, blue circles), merging clusters (n m /™, 
red squares), and scintillation clusters (n S ci/c C \ green trian- 
gles). The superposition of these components is also shown 
(purple inverted triangles). 

The total number of reconstructed electrons after ap- 
plying both of elD-Cut and npe-Cut for the three clus- 
ter types, which are represented as h a , h m and n sc i, 
are calculated by applying the npe-Cut efficiencies of 
C™dqf»(q), f££dqf?{q) and dqf s c "(q) to the 

fit results, n s , n m and n SC j, respectively. In the integrals, 
q m - m and g max represent the HBD charge boundaries in 
the npe-Cut of 8 p.e. and 28 p.e. (4 p.e. and 17 p.e. for 
the low-gain sector). The variables, n, are also summa- 
rized above. Figure [5] shows the yield spectra from the 
calculation as functions of pt- 



2. Yield estimation of separated photonic electrons 

The estimated n s is the sum of nonphotonic electrons 
and photonic electrons which create the separated clus- 
ters in the HBD. In the following description, we denote 
the photonic electrons which create merging clusters as 
merging photonic electrons (MPE) and those which cre- 
ate separated single clusters as separated photonic elec- 
trons (SPE). In this section, the number of SPE is esti- 
mated to obtain the yield of the nonphotonic electrons. 




FIG. 7: (color online). A half of an annular region around the 
reconstructed electron track on the HBD for the definition of 
hbdringcharge. The inner and outer radii of the annular 
region are 7.0 cm and 8.0 cm respectively. The direction of 
the half region is determined as the opposite side to the edge 
of the HBD sector to avoid inefficiency around the edge. 
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FIG. 8: (color online). HBD charge distribution in the an- 
nular region for the reconstructed electrons with a transverse 
momentum ranging from 0.75 GeV/c to 1.00 GeV/c (solid 
black line), and the fitting result of the charge distribution 
for electrons with correlated charges (n spe /r PC , green circle) 
and without correlated charges (n non _ a p / r non ~ spc , blue square), 
and the superposition of the fitting results (red inverted tri- 
angle) . 



In the case where a reconstructed electron track is iden- 
tified as an SPE, the partner electron generates an addi- 
tional signal in the HBD, as illustrated in Fig. 4(c)| This 
property is utilized to estimate the number of SPE. For 
this estimation, we defined a new value, g r ing, as 



hbdringcharge: <7 r i ng The total charge in the HBD 
pads centered on a half of an annular region with an 
inner radius of 7.0 cm and an outer radius of 8.0 cm 
around the track projection of HBD as shown in 
Fig. [7] To avoid inefficient regions around the edges 
of the HBD sectors, we use one half of an annular 
region oriented away from the nearest sector edge 
(see Fig. E]). The 

9ring value is normalized by the 
area of the half of the annular region in the defini- 
tion. 

The choice of 7.0 cm to 8.0 cm is determined by three 
factors: (1) the distribution of distance between sepa- 
rated clusters of SPE has a maximum around 7.0 cm, 
(2) few HBD clusters have radii larger than 7.0 cm, 
and (3) larger area includes more scintillation back- 
ground and decreases the signal to background ratio. 
Whereas the (fring distributions for the nonphotonic elec- 
trons and MPE comprise signals only from scintillation 
light, the distributions for SPE include the correlated sig- 
nals around the tracks in addition to scintillation light. 

The variables used in the hbdringcharge analysis are: 

jspe p ro bability distribution of (fring for SPE 

jnon-spo p ro bability distribution of q r i n g for nonphotonic 

electrons and MPE 
n spo Number of SPE after applying the elD-Cut and npe- 

Cut. 



"■non-spo Number of electrons other than SPE after ap- 
plying the elD-Cut and npe-Cut. 

The /^(fcng) for SPE and the /? ww * ((fa n g) for non- 
photonic electrons and MPE can be estimated by hadron 
tracks and electron tracks with large q c \ us values, which 
consist almost entirely of MPE. Because hadrons and 
MPE clusters do not create any correlated signals around 
their tracks, the q r i ng distributions of the tracks are cre- 
ated by only the scintillation light. 

The / r spe (g r in g ) was estimated by using simulations. 
The dominant photonic electrons come from the Dalitz 
decays of tt° and rj and 7 from their decays which con- 
vert in materials. We simulated the detector responses 
for the Dalitz decay and the 7 conversion events of the 
neutral mesons by a GEANT3 simulation [35[ configured 
for the PHENIX detector system. The ir° and 77 spec- 
tra were parametrized in the simulation by mr-scaled 
Tsallis distributions (36j . together with their known 
branching ratios to Dalitz decays and 7 decays. In 
order to include contributions from scintillation light, 
/ r non_spc faring), which is identical to the Qring distribution 
from only the scintillation light, was convoluted to the 
result to obtain / r spo faring)- 

The 

9ring distribution for the reconstructed electrons 
selected by applying elD-Cut and npe-Cut was fit- 
ted with the superposition of the <? r ing distributions, 
/ r spe fa rin g) and /r n ~ sp °(<7ring), as 



"-spc X / r Spe (g rin g) + Tl 



non-spc 



x/ t n 



faring), (5) 



wher 



^spc 



and 



^non-spc 



are fitting parameters and repre- 



sent the numbers of SPE and other electrons in the q r i ng 
distribution, respectively, as summarized above. Similar 
to the <7cius distribution, the fitting for the (? r ing distribu- 
tion was also performed for each electron px region and 
each HBD sector. Figure [8] shows a fitting result in one 
HBD sector in the electron px region from 0.75 GeV/c 
to 1.00 GeV/c. 



3. Yield estimation of heavy flavor electrons 

Using the above fitting results of h s and n spe , the yield 
of nonphotonic electrons, iV npc was estimated with the 
formula 



N Dpc { Pl 



n s {p T ) 



Apt)- 



(6) 



The remaining background for the heavy flavor electrons 
in the nonphotonic electron sample comes from K e % de- 
cays and e + e~ decays of light vector mesons, namely p, 
u>, and <f>. Electrons from the Drell-Yan process also con- 
tribute to the background, however the contribution is 
known to be less than 0.5% of total heavy flavor electrons 
in this pt range and can be ignored. We determined the 
yield of the heavy flavor electrons from N npe by subtract- 
ing the components of the K e % electrons, which are es- 
timated by simulation using a measured K cross section 
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TABLE I: Relative systematic uncertainties given on the 
heavy flavor electron yield. 



source 


uncertainty 


Pt range 


(GeV/c) 


hbdringcharge fitting 


16% 


( 0.50 < p T 


< 0.75 ) 




6% ~ 4% 


( 0.75 < pt 


< 1.75 ) 




2% 


( 1.75 < p T 


) 


hbdcharge fitting 


2% 


( 0.50 < p T 


< 0.75 ) 




< 1% 


( 0.75 < p T 


) 




4% 


( 0.50 < p T 


< 0.75 ) 




< 1% 


( 0.75 < p T 


) 


hadron misID 


4% 


( 0.50 < p T 


< 0.75 ) 




< 1% 


( 0.75 < pt 


) 



|36| . and the electrons from light vector mesons, which 
are already estimated in previously published result [24| , 
as 

N HFc ( PT ) = N^ C ( PT ) - N Ke3 ( PT ) - N LVM (p T ), (7) 

where N K£3 (pt) and N lvm (pt) represent the electrons 
from the K e ^ decays and the light vector meson decays 
respectively. 



D. Systematic Uncertainty 

The systematic uncertainties for the heavy flavor elec- 
tron yield come from the fits for the <7 c ius distribution 
and the feng distribution, and from estimations of K e3 
contribution and misidentified hadrons. 

The most significant source in these contributions is 
the fitting uncertainty for the (fring distribution. We var- 
ied the radius of the annular region to an inner radius 
of 6.0 cm and an outer radius of 7.0 cm and also to 
8.0 cm and 9.0 cm from the default radii of 7.0 cm and 
8.0 cm. The uncertainty from the fitting was set to the 
amount of variation in 



after these changes. The 
estimated uncertainties decrease from about 16% of the 
heavy flavor electron yield in the momentum range of 
0.50 <p T < 1.00 GeV/c to about 2% above 1.75 GeV/c. 

The fitting uncertainty for the q c ius distribution comes 
from the estimation of the bias in the charge distribu- 
tion shape due to the electron pair reconstruction. The 
systematic uncertainty from this effect is estimated to be 
less than 2% by simulation. 

In the low momentum region, 0.50 < pr < 1-00 GeV/c, 
uncertainties from the K e % contribution and the hadron 
misreconstruction are significant. The uncertainty from 
the -fT e 3 contribution comes almost entirely from the un- 
certainty on the K cross section used in the K e ^ simula- 
tion. This uncertainty amounts to about 4% of the total 
heavy flavor electron yield in the low momentum region 
and decreases to less than 1% for px > 0.75 GeV/c. We 
also estimated the upper limits of the hadron contam- 
ination due to misreconstructions employing a hadron- 



enhanced event set. As a result, we determined the up- 
per limits as 4% of the total heavy flavor electron yield in 
the low momentum region which decreases to less than 
1% over 1.5 GeV/c. The upper limits are assigned as 
the systematic uncertainties from hadron misreconstruc- 
tions. Table U summarizes the systematic uncertainties 
on the heavy flavor electron yield. 



E. 



Results of Heavy Flavor Electron Yield and 
Signal Purity 



From Eq. [5] and Eq. [7] and the discussion in Sec. IHI Dl 
the heavy flavor electron yield spectrum with the sys- 
tematic uncertainties was determined. The spectrum is 
shown in Fig. [9] We also show the yield of inclusive 
reconstructed electrons after applying the elD-Cut and 
npe-Cut and the estimated K e % contribution. The elec- 
trons from e + e~ decays of the light vector mesons are 
not shown in Fig. [9j but they are less than 5% of the 
heavy flavor electron yield in this px range. 
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FIG. 9: (color online). Heavy flavor electron yield spectrum. 
The black square points represent the total number of the 
reconstructed electrons after applying the elD-Cut and npe- 
Cut. The red circle points represent the estimated yield of 
the heavy flavor electrons. The yellow bands represent the 
systematic uncertainties for the heavy flavor electron yield. 
The green triangle points with dashed lines represent the esti- 
mated K e 3 contribution with systematic uncertainties shown 
by light-blue bands. 

The ratio of the nonphotonic electron yield to the pho- 
tonic electron yield in this measurement, 



R(pr) = 



N n P e (p T ) 



N* eco (p T ) - N n P c (p T ) 



(8) 



where N™ 00 denotes the total number of reconstructed 
electrons after applying the elD-Cut and npe-Cut, is 
shown in the top panel of Fig. [TUJ In Eq. [FJ we assumed 
the fraction of misidentified hadrons in the reconstructed 
electrons after the cuts is negligible as shown in Fig. UJ 
and so the number of photonic electrons can be repre- 
sented as Nl ec °(pT) — N upc (pt)- The same ratio from 
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TABLE II: Relative systematic uncertainties on the cross sec- 
tion due to uncertainties in the total sampled luminosity, 
trigger efficiencies, and detector acceptance. These system- 
atic uncertainties are globally correlated in all pr regions 
(pt > 1.25 GeV/c for the uncertainties on rig B ) - 

source uncertainty 
MB trig, cross sect. 9.6% 
acceptance A 8% 
reco. efficiency e rcc 6% 
MB trig, efficiency eJ£P 2.5% 
e trig, efficiency 3.6% 



in p T > 1.25 GeV/c 



a previous measurement [24| is also shown in the figure. 
The previous measurement employed two other methods 
for the background estimation, namely a cocktail method 
and a converter method. In the cocktail method, a sum 
of electron spectra from various background sources was 
calculated using a Monte Carlo hadron decay generated. 
This sum was subtracted from the inclusive electron sam- 
ple to isolate the heavy flavor contribution. With the 
converter method, a photon converter around the beam 
pipe was introduced to increase the photon conversion 
probability by a well-defined amount, and thus allow de- 
termination of the photonic background. The nonpho- 
tonic to photonic electron ratio is improved by a factor 
of about 2 or more in pt > 1.0 GeV/c compared with 
the previously measured result due to the rejection of 
photonic electrons by the HBD. 

The signal purity is defined as the ratio of the yield of 
the heavy flavor electrons to the reconstructed electrons 
after applying the elD-Cut and npe-Cut, 



D( PT ) = 



_ N (Pt) 



N™ co (pr)' 



(9) 



The result is shown as the bottom plot in Fig. [lOj We 
also show the result of the signal purity in the previous 
measurement. Comparing with the previously measured 
result, the signal purity is improved by a factor of about 
1.5 in a p T range from 0.75 GeV/c to 2.00 GeV/c. 



IV. HEAVY FLAVOR ELECTRON CROSS 
SECTION 

The invariant cross section is calculated from 



E- 



d 3 <i 



1 1 



N(A PT , Ay) 



dp 3 2np T L 74e rcc e t ri g Ap T Ay 



(10) 



where L denotes the integrated luminosity, A the accep- 
tance, e re c the reconstruction efficiency, etrig the trigger 
efficiency, and N the estimated number of heavy flavor 
electrons. 

The luminosity, L, was calculated from the number of 
MB events divided by the cross section for the MB trig- 
ger. For the latter, a value of 23.0 mb with a systematic 
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FIG. 10: (color online), (top) Ratio between the yields of 
the nonphotonic electrons and the photonic electrons in the 
reconstructed tracks. The red circles and the blue squares and 
triangles represent this analysis result and the previous result 
[24| . respectively. The error bars and bands represent the 
statistic and the systematic uncertainties, (bottom) Signal 
purity which is a ratio of the yield of the heavy flavor electrons 
to the total reconstructed electrons. 



uncertainty of 9.6% was estimated from van-der-Merr 
scan results [37J corrected for the relative changes in the 
BBC performance. The combination of the acceptance 
and the reconstruction efficiency, Ae roc (pT), was esti- 
mated by a GEANT3 simulation. We found that Ae rec (pT) 
has a value of 4.7% x (1±8 x 10- 2 (acc.)±6 x 10~ 2 (reco.)), 
with a slight px dependence. 

The efficiency of the MB trigger for the hard scattering 
processes, including heavy flavor electron production, is 
e{>g = 79.5% x (1 ± 2.5 x 10" 2 ). The efficiency of the 
electron trigger for the electrons under the condition of 



the MB trigger firing, e^ B (p T ) 



etrig (pt Ver- 



mes 



be calculated by the ratio of the number of the recon- 
structed electrons in the MB triggered sample in coin- 
cidence with the electron trigger to the number of the 
reconstructed electrons without the coincidence. The ef- 
ficiency is shown in Fig. QT] as a function of pt- 
Whereas we used the calculated efficiency values for the 
momentum region of px < 1-25 GeV/c, we assumed a 
saturated efficiency for pt > 1.25 GeV/c and estimated 
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the value with a fitting as shown in Fig. [TT] The fit- 
ting result is e p i a toau = 56.5% x (1 ± 3.6 x 10~ 2 ). The 
total trigger efficiency etrig(PT) can be calculated with 

the above two efficiencies as etrig(pr) = el-rig x e trig B (? , r)- 
Table [TT] summarizes the systematic uncertainties on the 
cross section due to uncertainties in the total sampled 
luminosity, trigger efficiencies, and detector acceptance. 
All systematic uncertainties listed in Table [TTJ are glob- 
ally correlated over whole pr region (pr > 1.25 GeV/c 



for the uncertainties on e. 



1MB 



)• 



0.8 
0.7 
0.6 



e2 0.5 



0.4 
0.3 
0.2 
0.1 





t/ndf= 7.04/7 
Plateau Eff. = 56.5 %x(]± 3.6xl0 2 ) 



0.5 



1.5 



2.5 



3.5 



4 4.5 5 
P T [GeV/c] 



FIG. 11: (color online). Efficiency of the electron trigger 
for reconstructed electrons under the condition that the MB 
trigger was issued. The red line represents the fitting result 
with the constant function and the green band represents the 
fitting uncertainty. 



The measured cross section of heavy flavor electrons 
is shown in Fig. [12] and tabulated in Table Mil A cor- 
rection for bin width |38| is applied to the pt value of 
each point. The figure also shows the previously pub- 
lished result (24|. The new result agrees well with the 
previous result within the uncertainties. Note that in 
this paper we employed a new analysis method with the 
HBD whereas the previous measurement employed dif- 
ferent methods, the cocktail method and the converter 
method. The consistency between these measurements 
proves that additional photonic backgrounds generated 
in the HBD material are removed, and that this new 
analysis method with the HBD is robust. 

The electron cross section from J/ijj — > e + + e~ de- 
cays estimated by the cocktail method [25| and a fixed 
order next-to-leading log (FONLL) pQCD calculation of 
the heavy flavor contributions to the electron spectrum 
[39j are also shown in Fig. [12] The J/tp contribution 
to the heavy flavor electrons are less than 2% in pr < 
1.25 GeV/c and increase to -20% until p T = 5.0 GeV/c. 
The FONLL pQCD calculation shows that the heavy fla- 
vor electrons in the low momentum region are dominated 
by charm quark decays, and the contribution from bot- 
tom quarks in pt < 1.25 GeV/c is less than 5%. 
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FIG. 12: (color online), (top) Invariant differential cross sec- 
tions of electrons from heavy- flavor decays. The red circles 
are this analysis of 2009 data and the blue squares are the 
previous 2005 data [24|] for the nonphotonic electron cross 
sections. The error bars and bands represent the statistical 
and systematic uncertainties. The scaling uncertainty from 
the Vernier scan is not included in the error bands because 
the same uncertainty must be considered for both the results 
of 2009 and 2005. The purple dashed dotted line is electron 
cross section from J/iji — » e + + e~ decays estimated from 
the cocktail method 25] . The solid and dashed curves are 
the FONLL calculations, (bottom) Difference of the ratio of 
the data and the FONLL calculation from 1. The upper and 
lower curve shows the theoretical upper and lower limit of the 
FONLL calculation. 



V. HEAVY FLAVOR ELECTRON SPIN 
ASYMMETRY 

Since parity is conserved in QCD processes, thereby 
disallowing finite single spin asymmetries, using Eq. [3] 
we express the expected electron yields for each beam- 
helicity combination as 



N?I(N ,A LL ) = AT (1 + \P B Py\A LL ) 



r oxp/ 



Nl^(N ,A LL ) 
N?I(N ,A LL ) 
N e _^(N 0l A L 



ILL 



N (l + \P B Py\A LL )/r— 
N (l - \P B PY\A LL )/r + _ 
N (l-\P B P Y \A LL )/r_ + , 



(11) 



where N±±(Nq,All) denote the expected yields for col- 
lisions between the blue beam-helicity (±) and the yel- 
low beam-helicity (±) and Nq is the expected yield in 
collisions of unpolarized beams under the same inte- 
grated luminosity as the H — h beam-helicity combination. 
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TABLE III: Data table for the cross section result correspond- 
ing to Fig. m 



in the beam-helicity combinations, 



I-' i- 

V / C| 


dp^ 




stat. error 
[mbxGeV- 2 c 3 ] 


syst. error 


U.OIZ 


Z. 1Z X 1U 


3 


0.04x10" 


-3 


0.47x10" 


3 


U.O04 


i .yo x iu 


-4 


0.09x10" 


-4 


1.11x10" 


■4 


i 1 1 p; 

1. 110 


Z./oX 1U 


-4 


0.03x10" 


-4 


0.37x10" 


-4 


l.ODD 


i no v in - 
i.uy x iu 


-4 


0.02x10" 


-4 


0.13x10" 


-4 


l.Dl / 


4. / / X IU 


-5 


0.08x10" 


-5 


0.58x10" 


-5 


1 Qfil 
l.OD / 


n Q/1 w 1 O" 
z.o4X IU 


-5 


0.05x10" 


-5 


0.27x10" 


-5 


O 1 1Q 

z. llo 


1. ID X IU 


-5 


0.04x10" 


-5 


0.13x10" 


-5 


z.ooy 


D.UoX IU 


-6 


0.20x10" 


-6 


0.68x10" 


-6 


o £1 o 
z.Dly 


o.zo X IU 


-6 


0.19x10" 


-6 


0.37x10" 


-6 


z.ooy 


i qo v in~ 
l.oz X iu 


-6 


0.11x10" 


-(j 


0.20x 10" 


-6 


q i on 


l.Uo X IU 


-6 


0.07x10" 


-6 


0.12x10" 


-6 


3.370 


6.20x10" 


-7 


0.41x10" 


-7 


0.69x10" 


-7 


3.620 


4.07x10" 


-7 


0.26x10" 


-7 


0.45x10" 


-7 


3.870 


2.42x10" 


-7 


0.19x10" 


-7 


0.27x10" 


-7 


4.121 


1.59x10" 


-7 


0.15x10" 


-7 


0.18x10" 


-7 


4.371 


1.07x10" 


-7 


0.11x10" 


-7 


0.12x10" 


-7 


4.621 


8.02x10" 


-8 


1.11x10" 


-8 


0.89x10" 


-8 


4.871 


5.38x10" 


-8 


0.71x10" 


-8 


0.60x10" 


-8 



TABLE IV: Systematic uncertainties by type. The scaling un- 
certainty denotes an uncertainty on scaling of the raw asym- 
metry A^ BG and the offset uncertainty denotes an uncer- 
tainty on the absolute value of the asymmetry. The "global" 
in this table means the uncertainties are globally correlated in 
all pr regions. The scaling uncertainty is represented as the 
ratio of the uncertainty to the signal (AS/S) given in percent 
and the offset uncertainty is represented as the absolute value 
of the uncertainty. 



uncertainty 



type 



signal purity D 
polarization ( — p-^rp ) 
relative luminosity r 
background asymmetry Af^ 



-6% 
8.8% 
0.14 x 10" 2 
0.2 x 10 



2 l-D 
x ~D 



scaling 
global scaling 
global offset 
offset 



N±±(Nq, All) are used for fitting functions to estimate 
All as described below. Pb and Py represent the polar- 
izations of the beams. The beam polarizations are mea- 
sured with a carbon target polarimeter [4fj| . normalized 
by the absolute polarization measured with a separate 
polarized atomic hydrogen jet polarimeter [4l|,[42| at an- 
other collision point in RHIC ring. The measured polar- 
izations are about P — 57% with a relative uncertainty 
of AP/P = 4.7 x 10~ 2 in the measurement. The relative 
luminosities are defined as the ratio of the luminosities 



r 

r+- 



= L ++ 

— L__ 

= L±t 
= L±± 



(12) 



where L±± represent the integrated luminosities in the 
beam-helicity combinations shown by the subscript. The 
relative luminosities are determined by the ratios of MB 
trigger counts in the four beam-helicity combinations. 

The double-spin asymmetry for inclusive electrons af- 
ter applying elD-Cut and npe-Cut, which include not 
only the heavy flavor electrons (S) but also the back- 
ground electrons (BG), is determined by simultane- 
ously fitting the yields of electrons in each of the four 
beam-helicity combinations with the expected values 
N±1(N ,A L l) from Eq.[[TJ where All and N are free 
parameters. To perform the fit, a log likelihood method 
assuming Poisson distributions with expected values of 
N*^±(No, All) was employed. The fit was performed for 
electron yields in each fill to obtain the fill-by-fill double- 
spin asymmetry. We confirmed that all asymmetries in 
different fills are consistent with each other within their 
statistical uncertainties and, therefore, the patterns of 
the crossing helicities in the fills do not affect the asym- 
metry measurement. The final double-spin asymmetry 
for inclusive electrons, A l + l bg {pt), was calculated as the 
weighted mean of the fill-by-fill asymmetries. 

The double-spin asymmetry in the heavy flavor elec- 
tron production, A^£ e , was determined from 



/t HFc 

A LL 



(Pt) 



1 



D( Pt ) 



aS+BG 
LL 



(Pt) 



1 - D(pt) 
D(pr) 



aBG 
LL 



(P'l 



(13) 

where A B( £ represents the spin asymmetries for the back- 
ground electron production, and D represents the sig- 
nal purity defined in Eq. [9] and shown in Fig. 1101 As 
previously discussed, most of the background electrons 
come from Dalitz decays of the ir° and 77, or from con- 
versions of photons from decays of those hadrons. The 
fractional contribution on the partonic level, and there- 
fore the production mechanism for the ir° and 
pected to be very similar up to ~ 10 GeV/c 
We assume identical spectra for double-spin asymme- 
tries of 7T° production and 77 production, and estimated 
A®2 from only the ir° double-spin asymmetry using data 
from this PHENIX measurement. The resulting A® 1 is 

2 in 0.5 < p T < 2.5 GeV/c 
2 in 2.5 < p T < 3.0 
GeV/c, with uncertainties less than 0.2 x 10~ 2 . 

Systematic uncertainties on A^£ e are separated into 
scaling uncertainties and offset uncertainties. The scal- 
ing uncertainties come from uncertainty in the beam po- 
larizations, Pb and Py, and the signal purity, D. The 
uncertainty from the beam polarization is estimated as 
A{P b Py)/P B Py = 8.8% which is globally correlated 
over the whole pt range. The offset uncertainties come 
from uncertainties in the relative luminosity, r, and the 



-0.1 xl0~ 2 < Af£ < 0.1x10" 
and 0.1 x 10~ 2 < Af2 < 0.2 x 10 
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background asymmetry, The uncertainty from rel- 

ative luminosity which is also globally correlated over is 
determined as Ar = 1.4 x 10~ 3 from comparison of the 
measured relative luminosities with the MB trigger and 
the ZDC trigger. The systematic uncertainties are sum- 
marized in Table HVl 

A transverse double-spin asymmetry Att, which is de- 
fined by the same formula as Eq. [5] for the transverse po- 
larizations, can contribute to All through the residual 
transverse components of the beam polarizations. The 
product of the transverse components of the beam po- 
larization is measured to be ~ 10~ 2 in this experiment. 
For 7T° production, the Att is expected to be ~ 1CT 4 
based on an NLO QCD calculation [43J. If we assume the 
transverse asymmetries of ir° and heavy flavor electrons 
are comparable, we arrive at the value of All ~ 10~ 6 . 
This value is negligible compared with the precision of 
the A^ J [ BG measurement of ~ 10 -3 . 

The result of the double-spin asymmetry of heavy fla- 
vor electrons is shown in Fig. Q2] and tabulated in Ta- 
ble [V] We show systematic uncertainties for scaling and 
offset separately in the figure. The measured asymmetry 
is consistent with zero. 



,0.15 



0.05 



-0.1 



-0.15 



A"[ e p+p\fs=200GeV flu |<0.35) 



T 



I syst. error from dilution factor 
syst. error from rel. him. and A BI 

(not including 8.8% pol. scale error) 



0.5 



1 



1.5 



2.5 



3 

P T [GeV/c] 



FIG. 13: (color online). Double-spin asymmetry of the heavy 
flavor electron production. The red error bands represent 
scaling systematic uncertainties from the dilution factor and 
the blue error bands represents offset systematic uncertainties 
from relative luminosity and the background spin asymmetry. 



VI. DISCUSSION 

In this section, we discuss constraint of Ag from the 
measured double-spin asymmetry with an LO pQCD cal- 
culation. In p+p collisions at y/s = 200 GeV, heavy 
flavor electrons with momentum ranging 0.50 < pt < 
1.25 GeV/c are mainly produced by open charm events 
as described in Sec. [IVJ Whereas the precise mechanism 
for J/ij) production is unknown, unpolarized and polar- 
ized cross section of the open charm production can be 
estimated with pQCD calculations. In LO pQCD cal- 
culations, only gg — > cc and qq — ¥ cc are allowed for 



the open charm production. The charm quarks are pri- 
marily created by the gg interaction in the unpolarized 
hard scattering. In addition, the anti-quark polarizations 
are known to be small from semi-inclusive DIS measure- 
ments precisely enough that both DSSV [19( and GRSV 
44j expect contribution of polarized qq cross section to 
the double-spin asymmetry of the heavy flavor electrons 
in \r)\ < 0.35 and p T < 3.0 GeV/c to be - 10~ 4 [Hj], 
which is much smaller than the accuracy of this mea- 
surement. Therefore, in this analysis of Ag, we ignore 
the qq interaction and assume the asymmetries are due 
only to the gg interaction. Under the assumption, the 
spin asymmetry of the heavy flavor electrons is expected 
to be approximately proportional to the square of polar- 
ized gluon distribution normalized by unpolarized distri- 
bution, \Ag/g(x,n)\ 2 . 



We estimated the unpolarized and the polarized cross 
section of charm production in p+p collisions with a 
LO pQCD calculation of gg — > cc 20J. For this cal- 
culation, CTEQ6M 45] was employed for the unpolar- 
ized parton distribution functions (PDF). For the polar- 
ized PDF, we assumed |Ag(x,/i)| = Cg{x,ii) where C 
is a constant. The charm quark mass was assumed as 
m c =1.4 GeV/c 2 and the factorization scale in CTEQ6 
and the renormalization scale were assumed to be iden- 
tical to \x = m T ': 



p T 2 + m c 2 . 



The fragmentation and decay processes were simulated 
with pythia8 [Ifllilj]. We generated pp — > cc + X events 
and selected electrons from the charmed hadrons, D + , 
D°, D s , A c and their antiparticles. We scaled the charm 
quark yield in pythia with respect to the pQCD calcu- 
lated unpolarized and polarized cross sections to obtain 
unpolarized and polarized electron yields from charmed 
hadron decays under these cross sections. We also ap- 
plied a pseudorapidity cut of \ij\ < 0.35 for the electrons 
to match the acceptance of the PHENIX central arms. 
The shape of the expected spin asymmetry A l 1 l ?c (pt) is 
then determined from the simulated electron yields. 



Figure [T4l shows the distributions of the gluon Bjorken 
x contributing to heavy flavor electron production in 
the momentum range 0.50 < pt < 1-25 GeV/c, from 
pythia. Using the mean and the RMS of the distri- 
bution for 0.50 < pt < 1-25 GeV/c, we determine 
the mean x for heavy flavor electron production to be 



(lo f 



10 • 



■1.6 



+0.5 



We calculat ed exp ected A l ^{pt) by varying C = 
\Ag/g\. Figure 15(a) shows several of these curves, along 
with the measured points, x 2 values are calculated for 
each value of C, along with related uncertainties. By as- 
suming that the systematic uncertainties on the points 
are correlated and represent global shifts, we defined the 
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TABLE V: Data table for the A™ e result corresponding to Fig. [H 



MGeV/c] 


Aft 


stat. uncertainty 


syst. uncertainty (offset) 


syst. uncertainty 


0.612 


2.83xl0~ 2 


2.66xl0 -2 


0.75xl0^ 2 


0.50X10" 2 


0.864 


-1.20xl0~ 2 


1.21xl0" 2 


0.30xl0~ 2 


0.08xl0~ 2 


1.115 


0.76xl0" 2 


1.30xl0~ 2 


0.21xl0" 2 


0.04xl0~ 2 


1.366 


2.08xl0~ 2 


1.63xl0~ 2 


0.18xl0" 2 


0.10xl0~ 2 


1.617 


-0.69xl0" 2 


2.18xl0 -2 


0.17xl0~ 2 


0.03xl0~ 2 


1.867 


-1.39xl0" 2 


2.68xl0~ 2 


0.16xl0~ 2 


0.03xl0" 2 


2.118 


4.82xl0" 2 


3.46xl0~ 2 


0.16xl0~ 2 


0.09xl0~ 2 


2.369 


-5.91xl0 -2 


4.40xl0 -2 


0.16xl0~ 2 


O.llxlO -2 


2.619 


-6.97xl0~ 2 


5.47xl0~ 2 


0.16xl0" 2 


0.13xl0" 2 


2.869 


6.43xl0~ 2 


7.07xl0~ 2 


0.16xl0~ 2 


0.12xl0~ 2 



-x10l_ 
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FIG. 14: (color online). Bjorken x distributions of gluons con- 
tributing the heavy flavor electron production with momen- 
tum ranging 0.50 < pr < 1.25 GeV/c obtained from pythia 
simulation. The distribution is normalized with respect to the 
number of total generated charmed hadrons. 



quantity x 2 as 



x 2 (C) 

P(C) 



-21og((27r) f P(C) 
/ dpdqN(p)N(q)x 



^.i scale 
/syst 



-(i+g7 a y:t caIo )/(^;g)) 



A(P B P Y ) 
PbPy 



where N(X) denotes normal probability distribution, i.e. 
N(X) = l/V27rcxp(-A" 2 /2), n is the number of the 
data points and equal to three, and for the i-th data 
point, Xi is the pr value, yt is the All value, and £g tat , 



offset 
c syst 



and 



scale 
syst 



represent the statistical, offset system- 



atic and scaling systematic uncertainties, respectively. 
f{pT',C) denotes the expected All{pt) for the parame- 
ter of C = \Ag/g\. A(PbPy) is an uncertainty for po- 
larization mentioned in Sec. |Vj If we set the systematic 



uncertainties, e 



offset 
syst 



and 7 ( 



i scale 
syst 



fined x 2 is consistent with the conventional \ 2 - 

The resulting x 2 curve is shown in Fig. |15(b)| plot- 
ted as a function of C 2 = \Ag/g\ 2 because the curva- 
ture becomes almost parabolic. The minimum of x 2 , 
Xmin' i s located at \Ag/g\ 2 = 0.0 which is the bound- 
ary of \Ag/g\ 2 . Ax 2 = X 2 ~ Xmin = 1 an d 9 were uti- 
lized to determine la and 3cr uncertainties. With these 
criteria, we found the constraints on the gluon polar- 
ization are \Ag/g({log w x) , /i)| 2 < 3.3 x 10~ 2 (lcr) and 
10.9 x 10~ 2 (3cr). The constraints are consistent with 
theoretical expectations for Ag/g(x,(x) at (log 10 x) = 

— 1.6±g-| and fx = 1.4 GeV which are 0.006 from 

DSSV, '~ 0.016 from GRSV(std) and - 0.019 from 
GRSV(val) using CTEQ6 for the unpolarizcd PDF. 

The effects of the charm quark mass and scale fac- 
tor in the cross section calculation were also checked 
by varying the charm mass from m c = 1.3 GeV/c 2 to 
1.5 GeV/c 2 and the scale to /i 2 = 0.75m^ 2 and 1.5m^ 2 . 
Figure 15(b)| also shows the resulting x 2 curves. Consid- 
ering the variation of the crossing position at Ax 2 = 1, 
the constraint including the uncertainties from the charm 
mass and the scale can be represented as \Ag/g\ 2 < 
(3.3t^(mass)+°;4(scale)) x 10- 2 (1<t). 

The integral of the CTEQ6 unpolarized PDF in the 



sensitive x 



to zero, the newly de- 



region of (log 10 x) = — 1.6_Q4 and [i = 

1.4 GeV is J 001 dxg(x, /i) = 4.9. Hence the constraint 
on the integral of the polarized PDF at la corresponds 
to | Joo] 8 dxAg(x, < 0.85. This study also highlights 
the possibility for constraining Ag in this Bjorken x re- 
gion more precisely in the future with higher statistics 
and higher beam polarizations. 



VII. SUMMARY 

We have presented a new analysis method for identi- 
fying heavy flavor electrons at PHENIX. With this new 
method, the signal purity is improved by a factor of about 
1.5 around 0.75 pr ~ 2.00 GeV/c due to the rejection 
of photonic electrons by the HBD. We have reported 
on the first measurement of the longitudinal double-spin 
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(a)Expectcd for several \Ag/g\. 
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IAg/gl' 2 
(b)x 2 curves. 

FIG. 15: (color online), (a) Aft for \Ag/g\ = 0.00, 0.10, 
0.20, 0.30, 0.40, 0.50 are shown as the solid line, the dashed 
line, the dotted line, the dashed dotted line, the long-dashed 
dotted line, the dashed triplicate-dotted line respectively. 
They are plotted with the measured data points and the no- 
tation for the error bars are same as Fig. 1131 (b) % 2 curves 
calculated from (a) as a function of \Ag/g\ . The black solid 
line is the default configuration. The blue curves are after 
changing the charm mass to 1.3 GeV/c 2 (dashed line) and to 
1.5 GeV/c 2 (doted line) and the red curves are after chang- 
ing the scale /i 2 to 0.75my 2 (dashed dotted line) and 1.5m^ 2 
(long-dashed dotted line). 

asymmetry of heavy flavor electrons, which are consistent 
with zero. Using this result, we estimate a constraint of 



\Ag/g(\og 10 x = -1.6±°j,/i = m c T )\ 2 < 3.3 x 1(T 2 (1ct). 
This value is consistent with the existing theoretical ex- 
pectations with GRSV and DSSV. With improved statis- 
tics and polarization, the hclicity asymmetry of heavy 
flavor electron production can provide more significant 
constraints on the gluon polarization, and complement 
other measurements of AG. 
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